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NASTRAN Modifications for 
Recovering Strains and Curvatures 

INTRODUCTION, 

The Lewis Modifications Project has two basic objectives: 

1. Produce strain and curvature values for the general two-dimensional 
elements (TRIAl, TRIA2» QUADl and QUAD2) as an alternative to the 
currently produced stress values. 

2. Transform the strain/ curvature (or stress) values to a "material'* 
coordinate system and interpolate to the grid points to which the 
elements are connected. 

For the purpose of this project, the term strain/curvature means the 
quantities defined by 
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where (u, v, w} is the displacement vector measured in a local Cartesian 
rectangular x, y, z coordinate system. 

This report describes the programming changes made to implement 
these objectives and provides usage instructions for the engineer using 
this new capability. 
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PROGRAMING CONSIDERATIONS 


All modifications to NASTRAN are upward compatible in the sense that no 
change will be observed by any user of currently defined Level IS. 5 capability. 

Part 1 

Part I consists of program modifications needed to produce strain/ 
curvature values in lieu of stress values for the general two-dimensional 
elements. This task was accomplished by a modification of routines STRBSl, 
STRMEl, and STRQD2. Mathematically, the stress- strain matrix was replaced 
by an identity matrix, and thermal effects were removed. Strains resulting 
from thermal expansion arc included in the displacement-related computations. 
The strains and curvatures are returned by the element routines in the data 
words that were used for stresses. Then the output quantities are 
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The strain and curvature computations occur only if DIAG 23 has been turned 
on by the user. Otherwise, the modified code has no effect and the usual 
stresses are output. The output formats were not changed; interpretation of 
the output is simply accomplished by direct correlation with the DMAP 
sequence utilized; therefore, no work was required in any other module. 

Part II 

Part II consists of program modifications needed to transform the 
strain/curvature (or stress) values for the elements to a common "material” 
coordinate system and then interpolate to the grid points to which the 
elements are connected. A number of minor modifications to existing modules 
and a new module (CURV) were required to accomplish this task. Tliese will be 
described in some detail in the sections that follow. 

IFP Modificati o ns 

The MATl and MAT2 cards were modified to allow the inclusion of a coord- 
inate system reference identification number. For the MATl card, MPT now 
has 12 entries, the twelfth one of which is a non-negative integer coordi- 
nate system reference identification number. For the MAT2 card, a new 
third continuation card holds the 17th data item, which means the same 
thing as the 12th entry on the MATl card. 

Material Subroutine Modifications 

The MPT data block contains the coordinate system reference identifica- 
tion number. No functional modification to the material property processing 
routine PREMAT was necessary. No change is necessary for any module which 
does not desire to use the new data item on the MATl and MAT cards. 


XMPLBD Modifications 


A new entry was added for the new modulo. 

Module Name: CURV 

Subroutine Name: CURV 

MPL Properties: 

Inputs - 6 

Outputs •* 2 

Scratch - 4 

Parameters - 3 

1 - Integer, default = 0 

2 - Integer, default = 0 

3 - Integer, default = 0 

CURV 

Anew functional module, CURV, produces strain/curvature (or stress) 
data at the element grid point locations. A DMAP Alter is used by the 
user at execution time to invoke this module. The MPL characteristics follow. 


Input Data Blocks 


1 

0ES1 

0FP Element strain/curvature (or stress) table 

2 

MPT 

Material Properties Table 

3 

CSTM 

Coordinate System Transformation Matrices 

4 

EST 

Element Summary Table 

5 

SIL 

Scalar Index List 

6 

GPL 

Grid Point List 

Output Data Blocks 


1 

0ES1M 

0FP Element Strain/Curvature (or stress) 
Table in Material Coordinates 

2 

0ES1G 

0FP Grid Point Strain/Curvature (or stress) 
Table 


Note: If DIAG 23 is on, strain /curvature values are computed; otherwise, 

stresses are generated. 
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Scratch Data Blocks 


1 SCRATCHl 

2 SCRATCH2 

3 SCRATCH3 

4 SCRATCH4 


used to hold various 
lists and tables 
during the execution 
of the module. 


Parameters 

1 Integer, Input, default » 0, Output option flag 

0 - pass thru the print/punch/plot device 
code from 0ES1 to 0ES1G 

- use value for 0ES1G output device code 

2 Integer, default * 0, Processing option flag 

0 - generate both 0ES1G and 0ES1M 

1 - generate only 0ES1M 

3 Integer, default = 0, Number of interoolation points 

0 - use all elements in the interpolation 

j^O - use value to determine and use closest 
elements in the interpolation. 


CURV Module Data Processing Operations 

Each subcase is processed independently of all other subcases, one at a 
time, until the 0ES1 data block is exhausted. For each subcase, two phases of 
processing occur as described below in summary and later in detail. Either 
stresses or strain/curvatures are computed according to the content of 0ES1. 

A logic flow chart of the CURV module is shown in Figure 1. 

CURV Phase 1 - Data collection and transformation to material coordinates 

Several tasks occur during Phase 1, siime of which only occur if P2=0. 

1. Build an abbreviated subset of the EST data block Ccalled ESTX) 
containing all elements of potential interest. The order will be 
that of the 0ES1 data block in order to simplify subsequent processing. 
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Elements not selected for output in case control for a particular 
subcase, while they will bo placed on ESTX, will not bo used for 
that subcase, h{*te that, as in the case of element stresses, 
strain/ curvatures will only bo produced for elements selected in 
case control, usually via STRESSsALL. 

2. Build a preliminary table of independent points consisting of the 
center locations of the elements contained on ESTX. The final table 
is generated after the coordinate system matrices are loaded 

into core. 

3. Build a preliminary table of dependent points consisting of the 
locations of all grid points to which elements on ESTX arc connected. 
Tl»e final table is generated later when coordinate system 
information is available and after external grid point identifications 
are computed. 

4. Build a list of MCSIDs which are present on ESTX. Write this list 
onto SCRATCHl for subsequent use. 

5. Transform the strain/curvature (or stress) data to tl.e material 
coordinate system. This data is written on data block 0ES1M 
for subsequent processing by module 0FP. During this task, the 
final lists of independent and dependent points are generated. 

CURV Phase 2 - Projection and interpolation 

For each material coordinate system of interest, the following tasks 
are carried out: 

1. Collect the element center locations (independent points). 

2. Collect the unique set of connection point locations (dependent 
points). All location data is expressed in basic coordinates. 

3. Convert all location coordinates to a "local” system and select the 
appropriate mapping projection. Scale the mapping variables con- 
sistently. 

4. Form the required entry lists and call the SSPLIN routine to produce 
the interpolation transformation matrix G as described in P.M. 

Section 3.4.85 for each dependent point. 
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5. Transform the strain/curvaturo (or stress) vectors from the element 
centers to the grid points. 

6. Recompute invariants. 

7. Generate the 0FP data block 0ES1G for subsequent printing and/or 
punching. The data is sorted on grid point ID. 

The geometry being processed in shown schematically in Figure 2. 
Detailed programming descriptions of these steps are presented in the 
remainder of this section. 

Phase 1 - Step 1 

The abbreviated EST (ESTX) consists of one or more records (one 
record per element type) of a single word containing the element type code 
followed by several groups of 11 (for TRIAi elements) or 14 (for QUADi 
elements) words as follows: 

Word Symbol Item 

1 ELTYPE Eiomtjnt Type code 

2 EID Element ID 

3 MCSID Material Coordinate System ID 

4-6 Xj,yj,Zj 

7—9 

10-12 

13-15 ^4»Y4»*4 

(for QUADS only) 

In order to generate this information, the following tasks are 
carried out: 


Basic location coordinates 
of connectivity points (3 
for triangles, 4 for QUADs) 


group of 11 or 
14 words repeats 
for each element 
of the specified 
type 
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a. Read the MPT data block, extracting the material Identification 
number (MID) and the material coordinate system identification number (MCSID) 
for each MATl and MAT2 material. Entries having no MCSID are ignored. This 
paired list is sorted on MID and written onto SCRATCHl for future use. 

b. The EST data is passed through core. Only element types of interest 
ore examined. Elements referencing materials having no MCSID are ignored. 

For elements of interest, find the MCSID from the pair list create*) in task (a). 


Phase 1 - Step 2 

The preliminary independent point list consists of six words for each 
element on ESTX as follows: 


Word 


Item 


1 

2 

3 

4 

5 

6 


Material coordinate system identification number 

Will be filled in with coordinate system type 
code later. 

Element identification number 


Center basic location coordinates 


Phase 1 - Step 3 

The preliminary dependent point list consists of 16 or 20 words for 
each .element on ESTX as follows: 



Word 


1 
2 

3 

4 

5 

6 

7 

8 

Phase 1 - Step 4 

During generation of ESTX« the entries of the MID •* MCSID pair list 
were flagged when referenced. At this point, the unreferenced items in this 
list are removed. 

Phase 1 - Step 5 

The 0ES1 and ESTX data blocks are now simultaneously read, element by 
element. For those entries appearing on both data blocks, the strain/ 
curvature (or stress) data items are transformed to the material coordinate 
system via subroutine TRANEM. The invariant quantities are recomputed* and 
the resulting data are written onto 0ES1M for subsequent processing by 0FP. 
As a preprocessing task, once the CSTM data has been read into core, the 
independent and dependent lists are converted to final forva by substituting 
the coordinate system type code for the second word of each entry. This 
task is only done if P2«0. As a po3t>processing task, items not on 0ES1M 
are removed from the point lists. 

*See Section 4.87.4.6 of the NASTRAN Programmer's Manual 


Item 


Material coordinate system identification number 

Will bo filled in with coordinate system type 
code later. 

Element identification number 
3 for TRIAi elements, 4 for QUADi elements 
SIL for connected point 


Basic location 
coordinates 


repeats for 
each connected 
point 
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Tlio transformation matrix U was designed to transform "tensor” components 
such as stresses. If strains are to be transformed (DIAG 23 is on), then 
the third component (y) must first be multiplied by 4. The (o)^^ « U(a)^ 
material components are computed by matrix multiplication. The, if DIAG 23 
is on, the "engineering" shear strain is recovered by multiplying the last 
component by 2. 

Phase 2 - Step 1 (phase 2 is done only if P2 » 0) 

Pass the independent point table, eliminating any points not associated 
with MCSID. Since this table has already been reduced to 0ES1M size, the 
resulting subset is exactly the desired data required in Phase 2. 


Phase 2 - Step 2 


Pass the dependent point table, eliminating any points not associated 
w^.th MSCID. Collect the unique set of grid point ID's and their corres- 
ponding basic location coordinate for subsequent use in Phase 2. This 
list will be sorted on external grid point ID. 

Phase 2 - Step 3 


a. Convert both the independent and* dependent location coordinates to 
the local material coordinate system. 


*^local 


t'^{e-v} 


E - r measured in basic coordinate system 

V - reference vector obtained from the CSTM data block (words 3-5) 

T - transformation matrix obtained from the CSTM data block (words 
(6-14) stored by row. 


See P.M. Section 2.3-33 and 3.4-64 for details. 


b. Convert dependent and independent lo'^tion coordinates {r^j^} ■ (x,y,z) 
to mapping coordinates {p} = as follows. For rectangular 

coordinate systems, {p} ■ {tp}. For cylindrical coordinate systems, 


p^ - r - 1 ^ 
“1 

mar 

P2 


. A . r > 0 

® c 0 , r - 0 


>3 - z 


as shown on the following sketch. 


(in radians) 

-IT < 0 < + 7T 


2 



- 11 - 


Por spherical coordinate f,ystems, {p} ■ {r, 0, <}>} where, if we let 


I ■ 



P 


Pj ■ t -/x^+y^+z^ 


P2 

P3 


, . { tan‘^(a,z), r > 0 

(in ridians)® , r • 0 

- 4, - i J 2 

(in ridians)® * *- - 0 


0 < 6 S TT 
-TT < (|> < TT 


as shown in the sketch below: 


z 
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c. Select the projection based on the Independent points. To do this we 
want to find the napping coordinate which has the smallest range of 
values. Define a characteristic length as follows for angular 
mapping coordinates: 

cylindrical: 1 “ average r 
spherical: X = average r 


The projection is selected by finding which coordinate yields 


' \in* ^max " ^min* ^max ” ^min^ 

for rectangular coordinate systems, rain “ ^®min* 

2 - z . } for cylindrical coordinate systems, or min (r„„^ - r„.„, 

mzx min ^ ’ max min’ 


Xe - A0 , 

max min 


, £d) - Xd) . } for spherical coordinate systems. 

' ^max ^min ^ ^ 


The resulting mapping surfaces are illustrated (for regular mesh 
spacing) on the following sketches. 


d. The interpolation independent variables are now selected by discarding 
the selected projection coordinatfis and scaling both the remaining 
two mapping coordinates to have the range >1 to 4-1. 

e. The interpolation dependent variables are now arbitrarily reduced 
and scaled by the same selection and scaling as was used for the 
independent variables. 

The net result of step 2 is two sets of pair lists, 

xi, yi, i = 1, Ni 
and xd, yd, d = 1, Nd 

which will now be used in step 4. 
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Cylindrical Mapping Surfaces 
(unifonaly spaced meshes are shown for convenience) 


r ■ const 


6 ■ const 


z 





Phase 2 ~ Step 4 

This stop consists of preparing the input data for interpolatica, calling 
the interpolation routine SSPLIN, and handling the resulting G matrix prepara- 
tory to carrying out Step 5, If the third parameter is nonzeroj its value is 
used to determine the number of closest independent points for each dependent 
point; only these points are used in the SSPLIN interpolation. If the third 
parameter is zero, all independent points are used in the SSPLIN interpolation 
to obtain all dependent points. In both cases, the points used are the ones 
shown on the previous mapping surface sketches. 

Phase 2 - Step 5 

Read into core the strain/curvature (or stress) data from 0ES1M and 
prepare it in a form suitable for multiplication by the transformation matrix 
G created in Step 4. Core is assumed to be large enough to hold the entire 
operation. Call GMATTS to carry out the transformation. 

Phase 2 - Step 6 

Generate the invariants associated with the strain/curvature (or stress) 
quantities generated in Step 5 using equations of Section 4.87.4.6 of the 
.NASTRAN Programmer's Manual. 

Phase 2 - Step 7 

From the data obtained in Steps 5 and 6, generate the 0ES1G data block. 
Tl»is data is generated in external grid point sort. 
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CURV 



0ES1G 



7. I Generate 0FP tables 


More Subcases^ 


Figure 1. CURV Module Logic Flow Chart 
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B 
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Grid Point Locations (Dependent Points) 
Element Center Locations (Independent Points) 
' Material Coordinate System 


Figure 2 , Element Geometry 



Subroutine TRANEM 


Purpose: Computes a transformation matrix U for triangles and quadrilaterals 

which will convert strain/curvature (or stress) vectors measured in 
the element coordinate system to a material coordinate system pro- 
jected on the surface of the element. 



coordinates 

Calling Sequence: 




[U] 


ielement 

coordinates 


CALL TRANEM(MCSID,NG,R,IC0MP,U,RC) 

MCSID - Material coordinate system identification number » input. 

NG - 3 for triangles, 4 5or quadrilaterals, input. 

R - basic coordinate locations of connection points, length»3*NG, input. 

IC0MP - 1 if material x-axis is used, 2 if material y-axis is used, output. 

U - transformation matrix, length=9, row stored, output. 

RC - basic location coordinates of element center, length»3, output. 

Requirement: The calling subroutine must call PRETRS. 

Method: 

1. Compute the element normal: 

Triangles - n » Vj^ x V23 

= CV3 - Vj) X (Vj - Vj) 


Quadrilaterals - n = Vjj x V24 (definition) 

■ (''s - ''l> “ ■ '2’ 
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2 . Compute the center of the element: 

3. Call TRANS to cor;;pute the unit vectors of the material coordinate system 
MCSID at the center of the element. 

4. Select the reference axis: 

If the element normal is not normal to the material coordinate x>axis 
projection, use the x-axis and set IC0MPsl. Othcntfisc, use the material 
coordinate y-axis and set IC0MP>2. The criteria is arbitrarily chosen 

to be n' ' > 0.4 for use of th>fe y-axis rather than the x-axis, 

m 

5. Compute the sine and cosine of the angle between the selected material 
coordinate system axis and the element coordinate system x-axis. 

6. Generate the transformation matrix U: 

Ujj ■ cos^0 
Uj 2 = sin^0 

= -sin 0 cos 0 
U 2 J ® sin^0 
U 22 = cos"0 
U 22 = sin 0 cos 0 
Ujj = sin 29 
Uj 2 = -sin 20 
» cos 20 
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USAGE OF SDR2 AND CURV MODULES 


General 

lliere are two basic choices available to the user. Firsts stresses, 
strains, or both may be requested. This is controlled by the use of DIAG 23 
in modules SDR2 and CURV. Second, the quantities may be computed in three 
basic ways; 

(a) At element centers using an element coordinate system. 

fb) At element centers using a material coordinate system. 

(c) At grid point locations using a material coordinate system. 

All combinations of (a), (b) and (c) are possible. This is controlled 
by the introduction of the CURV module » the second parameter of the module, 
and by selection of 0FP inputs in the ALTERS. Item (a) is the usual SDR2 
output, (b) is generated during the first phase of CURV, and (c) from the 
second phase of CURV. Output is controlled by CASE CONTROL (e.g., STRESS 
(PUNCH) -request) and by the second parameter of CURV. 

Data Requirements 

The first feature is the ability to^request strain/curvature data in lieu 
of the usual stress data for all TRIAl, TOUa, QUADl and QUAD2 elements selected 

for output. This feature is selected by turning on DIAG 23 in the Executive 

Control Deck. No other user data is needed for this feature. 

If both stresses and strains are desired for the selected set of elements, 
the following alter packet may be used. DIAG 23 is not used in the Executive 
Control Deck. 
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ALTER 121 

PARAM // C,N,SSST / C,N,+23 $ 

SDR2 CASECC,CS1W,MPT,DIT,EQEXIN,SIL,GPn.EDT,BGPDT, , ,UGV,EST,/ 

,,,0ES1A,,/ C,N, STATICS $ 

0FP 0ES1A,,,», // $ 

PARAM // C,N,SSST / C,N,-23 $ 

For this application^ element stresses will be writtn on data block 0ES1 and 
element strain/curvatures will be written on data block 0ES1A as a result of 
the alter. 

The second feature provides the user with the ability to output the 
element strain/curvatures (or stresses) generated according to the preceding 
discussion for all selected TRIAlt TRIA2, QUADl and QUAD2 elements in a 
material coordinate system at the grid points to which the elements are 
connected using surface spline interpolation. To use this feature, the 
following DMAP alter packet must be employed in displacement rigid format 
ntuitber 1, static analysis: 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST,SIL,GPL/0ES1M,0ES1G/C,Y,0UT0PT»O/C,N,O/C,N,O $ 
0FP 0ES1G,,,,,// $ 

0ES1M, while not subsequently used in the DMAP, must be present as it is used 
as a scratch file. 

The element strain/curvature (or stress) values may also be output in 
the material coordinate system. In this case, the tiilter would be 

alter 121 

CURV 0ES1,MPT,CSTM,EST,,/0ES1M,/C,Y,0UT0PT=O/C,N,1/C,N,O $ 

0FP 0ES1M,,,,,// $ 
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Both element and grid point related outputs may be obtained for the first 
CURV alter shown by changing the 0FP call to 

0FP 0ES1M,0ES1G,,,,// $ 

DMAP Module Description 

The calling sequence for the CURV module is described below: 

DMAP Call: 

CURV 0ES1,MPT,CSTO,EST,SIL,GPL / 0ES1M,0ES1G / C,Y,0UT0PT=O / 

C,Y,0PT2=O / C,Y,0PT3=O $ 

Input Data blocks: As described by standard documentation 

Note: SIL and GPL may be purged if 0PT2 0 . 

Output Data Blocks: 0ES1M - Element strain/ curvatures Cor stresses) measured 

in material coordinates 

0ES1G - Grid point strain/curvatures (or stresses) 

Note: 0ES16 may be purged if 0PT2 0 . 

Parameters : 

0UT0PT - Integer, Input, default = 0 . 

0 = Pass through tne Print/Punch/Plot device code from 0ES1 
to 0ES1M and 0ES1G., 

^0 - Use value for 0ES1M and 0ES1G output device code. 

0PT2 “ Integer, Input, default = 0 . 

0 - Generate both 0ES1G and 0ES1M 
r<0 - Generate only 0ES1M 

0PT3 ~ Integer, Input, default = 0 • 

0 - Use all independent points in the interpolation. 

/O “ Use value as the number of "closest" points to use in the 
interpolation. 
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The possible output device codes are: 


1 = print 

4 = punch 

5 = print and punch 


Output 

The case control card STRESS is used to request 
of output that are possible. Interpretation is made 


all of the various kinds 
by examining the values 


which appear under the heading labeled "FIBRE DISTANCE" as indicated below: 

i’ 


Quantity 


1. Element stress in 
element coordinates 

2. Element stress in 
material coordinates 


3. Grid point stress in 
material coordinates 


Quantities in "FIBRE DISTANCE" printout jj 

l; 

jj 

a, b, fibre distances Zj^ and Z 2 j| 

}> 

i| 

jr 

a. material coordinate system id (real) || 

b. axis code (real) | 

1.0 = material x-direction chosen i 

2.0 “ material y-direction chosen 

a, material coordinate system id (real) | 

b. projection code + 10*N (real) | 

1.0 = material x-axis is normal to projection jj 

2.0 = material y-axis is normal to projection 

3.0 = material z-axis is normal to projection j 

(N is the number of interpolation points used) j 


4. Element strain/curvatures a, 
in element coordinates 

5. Element strain/curvatures a. 
in material coordinates b. 


6. Grid point strain/curva- a. 
tures in material coordi-' b. 
nates 


b, fibre distances z^ and Z 2 [j 

li 

material coordinate system id (real) ! 

axis code (real) 

1.0 = material x-direction chosen | 

2.0 ' material y-direction chosen | 

material coordinate system id (real) | 

projection code + 10*N (real) j 

1.0 = material x-axis is normal to projection 

2.0 = material y-axis is normal to projection 

3.0 = material z-axis is normal to projection 

(N is the number of interpolation points used) 


{ 
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APPENDICES 


In order to assist the user in using this new capability, a sample 
problem is presented in Appendix A. In addition. Appendix B presents DMAP 
ALTER packets which illustrate how the various options may be combined to 
produce a desired result. The new user is encouraged to study these 
appendices before attempting to utilize this new capability in his work. 

Wliile Rigid Format 1, Static Analysis, was used in both appendices, 
the new capability described in this report is not restricted by Rigid Format. 
Any Rigid Format (or DMAP sequence for that matter) which produces real, 

S0RT1 stress output can be altered to use the new capability. In many cases, 
only the alter numbers will differ from the material presented herein. 
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APPENDIX A 


SAMPLE PROBLEM 


The data deck listing and selected NASTRAN output pages for the seall 
saa^)le problea shown below are presented on the following pages. The 
notes are keyed to the card nuabers on the deck listing. 
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Notes 


Card Nos . 
15 - 19 

23 

24 - 25 
26 

27 - 28 
29 

30 - 31 
32 

33 - 34 
47 - so 
54 


56 

58 

100 

111 


Description 

Generate and output stresses in the material coordinate system 
(0ES1M) and interpolate to the grid points (0ES1G). 

Turn on DIAG 23 so that all requests for STRESS will be treated 
as if they were calling for strain/curvatures. 

Compute element strains (0ES1A) . 

Turn off DIAG 23. 

Output element strains. 

Turn on DIAG 23. 

Generate element strain/curvatures in the material coordinate 
system (0ES1AM) and interpolate to the grid points (0ES1AG) . 

Turn off DIAG 23. 

Output element strain/curvature data. 

Global output requests which ^pply to all subcases. 

Request for all "stresses." Due to the nature of the DMAP alters, 
all of the items listed below will be printed as a result of this 
card. 

1. element stresses in the individual element coordinate systems. 

2. element stresses in the material coordinate system 

3. grid point stresses in the material coordinate system. 

4. element strain/curvatures in the individual element coordi- 
nate systems . 

5. element strain/curvatures in the material coordinate system. 

6. grid point strain/curvatures in the material coordinate system. 

Request for "stresses" for the elements (and connecting grid 
points) in set 2 to be punched. See card 54 description. 

Request for "stresses" for the elements (and connecting grid 
points) in set 3 to be printed and punched. See card 54 description. 

Defines coordinate system 1003 to be the material coordinate 
system associated with material 771. 

Defines the number of "closest" independent points to be used 
in interpolating to each dependent point. Used in the CURV DMAP 
instruction (see cards 17 and 31). 
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Sample Problem Element Strain/Curvatures 
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Sample Problem 

Grid Point Strain/Curvatures in Material Coordinates 
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Appendix B 
DMAP ALTER PACKETS 


Several modes of operation of using this new capability are possible. 
Some of those are listed below and described separately in the sections which 
follow. Combinations of the various modes of operation ore possible. The 
data requirements are combinations of the ones presented and can b / deduced 


from the examples given in the text. 


Element /A1 
Output I 
In h2 
Element I 
Reference\^3 


Standard Stress Output 

Generate strain/curvatures in place of stresses 
Generate strain/curvatures in addition to stresses 


Element /B1 
Output 1 
In <B2 
Material j 
Referencelps 


Generate stresses in material coordinates 

Generate strain/curvatures ,in material coordinates 

Generate both strain/curva?;ures and stresses in material coordinates 


Grid /Cl 

Output I 
In lC2 

Material | 
Referencc^pS 


Generate stresses at grid points (measured in material coordinates) 

Generate strain/curvatures at grid points (measured in material 
coordinates) 

Generate both strain/curvatures and stresses at grid points 
(measured in material coordinates) 


Element 
Output 
In Both 
Element S* 
Material 
Referenc 
(combina- 
tion of 


D1 


m 


'D3 


A5B) 


Both Ele-/El 
ment § GriA 
Output in *2 
Material j 
Reference iE3 
(comb in a- \ 
tion of B^C) 


Generate stresses in both element and material coordinates 

Generate strain/curvatures in both element and material coordinates 

Generate both strain/curvatures and stresses in both element and 
material coordinates 


Generate stresses at both element and grid point locations 
(measured in material coordinates) 

Generate strain/curvatures at both element and grid point locations 
(measured in material coordinates) 

Generate both strain/curvatures and stresses at both element and 
grid point locations (measured in material coordinates) 


Detailed data deck requirements will now be given for the modes of 
operation just enumerated. The standard output will be called option A1 and 
need not be described in detail. i^RECEDJNG PAGE jigrQj 


A2. To produce strain/curvature output insteitd of stress output for elements; 
Executive Control Deck: Turn on DIAG 23. 

Case Control Deck: Select desired elements via STRESS»requost. 

Bulk Data Deck: None 

Remarks; Only QUADl, (^AD2, TRIAl and TRIA2 elements will generate strain/ 
curvature output. Other elements will generate the usual stress 
values. 


A3. To produce botli stress output and strain/ curvature output for elements: 

Executive Control Deck: 

ALTER 121 

PARAM // C,N,SSST/ C,N,+23 $ TURN DIAG 23 0N 

SDR2 CASECC,CSTM,MPT,DIT,EQEXIN,SIL,GPTT,EDT,BGPDT,,,UGV,EST, / 

,,,0ES1A,, / C,N, STATICS $ 

0FP 0ES1A,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAH //C,N,SSST/C,N,-23 $ TURN DIAG 23 0FF 

ENDALTER 

Cass Control Deck: Select desired elements via STRESSarequest. 

Bulk Data Deck: None 

Remarks: 1. Only QUADl, QUAD2, TRIAl and TRIA2 elements will generate 

strain/curvature output. Other elements will generate the usual 
stress values in both sets of output. 

2. The same elements will appear in both the stress and strain/ 
curvature output. 
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Bl. To produce stress output measured in material coordinates: 


Executive Control Deck: 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,, // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTW,EST,. / 0ES1M, / C,Y,0UT0PT=O / C,N,1 $ 

0FP 0ESlM,,,^v U V,N,CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 

Case Control Deck; Select desired elements via STRESS= request. 

Bulk Data Deck: 

1, The material coordinate system identification number must be coded on 
the MATl or MAT2 card for any element whose stress is to be measured in 
material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QUADl, QUAD2, TRIAl and TRIA2 elements will appear in the stress output. 

2. A default basic material coordinate system is not provided. For this 
case, set up a rectangular system coincident with the basic system. 
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B2. To produce strain/ curvature output measured in material coordinates: 


Executive Control Deck: 


1 . 


2 . 


DIAG 23 


ALTER 

120,120 

0FP 

0UGV1,0PG1,0QG1,0EF1,, // V,N,CARDN0 $ 

ALTER 

121 

CURV 

0ES1,MPT,CSTM,EST,, / 0ES1M, / C,Y,0UT0PT=O / C,N,1 $ 

0FP 

0ES1M,,,,, // V,N,CARDN0 $ 

SAVE 

CARDN0 $ 

ENDALTER 



Case Control Deck: Select desired elements via STRESS=requcst . 

Bulk Data Deck: 

1. The mater;! al coordinate system identification number must be coded on 
the MATl or MAT2 card for any element whose stress is to be measured in 
material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QUADl, QUAD2, TRIAl and 'nUA2 elements will appear in the strain/ 
curvature output. 

2, A default basic material coordinate system is not provided. For this 
case, set up a rectangular system coincident with the basic system. 


B3. To produce both stress and strain/curvature output measured in material 
coordinates . 

Executive Control Deck: 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,, // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST., / 0ES1M, / C,Y,0UT0PT=O / C,N,1 $ 

0FP 0ES1M,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C.N.SSST / C,N,+23 $ TURN DIAG 23 0N 

SDR2 CASECC,CS1W,MPT,DIT,EQEXIN,SIL,GPTT,EDT,BGPDT,,,UGV,EST, / 

,,,0ES1A,, / C,N, STATICS $ 

CURV 0ESIA,MPT,CS1>1,EST,, / 0ES1AM, / C,Y,0UT0PT=O / C,N,1 $ 

0FP 0ES1/4M,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C,N,SSST / C,N,-23 $ TURN DIAG 23 0FF 

ENDALTER 

Case Control Deck: Select desired elements via STRESS=request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose stress and strain/ 
curvature is to be measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 
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Remarks : 


1. Only QUADl, QUAD2, TRIM and TRIA2 elements will appear in the 
stress and strain/ curvature output. 

2. A default basic material coordinate system is not provided. For 
this case, set up a rectangular system coincident with the basic 
system. 

3. The same elements will appear in both the stress and strain/ 
curvature output. 
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Cl. To produce stress output at grid points {measured in material coordinates) 

Executive Control Deck: 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,, // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST,SIL,GPL / 0ES1M,0ES1G / C,Y,0UT0PT»O $ 

0FP 0ES1G,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 

Case Control Deck: Select desired elements via STRESS=request . 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded on 
the MATl or MAT2 card for any element whose stress is to be measured 
in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only grid points connecting QUADl, QUAD2, TRIAl and TRIA2 elements 
will appear in the 0ES1G output. 

2. A default basic material coordinate system is not provided. For 

this case, set up a rectangular system coincident with the basic system. 

3. A working core of approximately words (where Nj^ and N 2 are the 

number of cells per side of a rectangular mesh) must be available to 
carry out the interpolation to the grid points for such a mesh. Sub- 
division of the surface mesh into sub-meshes may be necessary for large 
problems. 
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C2. To produce strain/curvature output at grid points (measured in material 
coordinates) ; 

Executive Control Deck; 

1. DIAG 23 

2. ALTER 120,120 

0FP 0UGV1,0PG1,PQG1,0EF1,, // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST,SIL,GPL / 0ES1M,0ES1G / C,Y,0UT0PT=O $ 

0FP 0ES1G,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 

Case Control Deck: Select desired elements via STRESS^^request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded on 
the MATl or MAT2 card for any element whose strain/ curvature is to be 
measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only grid points connecting QUADl, QUAD2, TRIAl and TRIA2 elements 
will appear in the 0ES1G output. 

2. A default basic material coordinate system is not provided. For 

this case, set up a rectangular system coincident with the basic system. 

3. A working core of approximately words (where Nj and N 2 are the 

number of cells per side of a rectangular mesh) must be available to 
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carry out the interpolation to the grid points for such a mesh. Sub 
division of the surface mesh into sub-meshes may be necessary for 
large problems. 



C3. To produce both strain/curvature and stress output at grid points 
(measured in material coordinates) : 

Executive Control Deck: 

Al.TER 120,120 

0FP 0UGV1 ,0PG1 ,0QG1 ,0EF1 , , // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTO,EST,SIL,GPL / 0ES1M.0ES1G / C,Y,0UT0PT=O $ 

0FP 0ES1G,,,,, // V,N,CARDN0 $ 

SAVE CARON0 $ 

PARAM // C,N,SSST / C,N,+23 $ TURN DIAG 23 0N 

SDR2 CASECC,CS1W,MPT,DIT,EQEXIN,SIL,GPTT,EDT,BGPDT,,,UGV,EST, / 

,,,0ES1A,, / C,N, STATICS $ 

CURV 0ES1A,MPT,CSTW,EST,SIL,GPL / 0ES1AM.0ES1AG / C,Y,0UT0PT<*O $ 

0FP 0ES1AG,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C,N,SSST / C,N,-23 $ TURN DIAG 23 0FF 

ENDALTER 

Case Control Deck: Select desired elements via STRESS»request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose strain/curvature or 
stress is to be measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 
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Remarks : 


1. Only grid points connecting QUADl, QUAD2, TRIM and 111^2 elements 
will appear in the 0ES1G and 0ES1AG outputs. 

2. A default basic material coordinate system is not provided. For 
this case, set up a rectangular system coincident with the basic 
system. 

3. A working core of approximately words (where Nj and N2 are 

the number of cells per side of a rectangular mesh) must be available 
to carry out the interpolation to the grid points for such a mesh. 
Subdivision of the surface mesh into sub-meshes may be necessary for 
large problems. 
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Dl. To produce stress output in both element and material coordinates: 


Executive Control Deck: 

ALTER 121 

CURV 0ESI,MPT,CSTW,EST,, / 0ES1M, / C,Y,0UT0PT=O / C,N,1 $ 

0FP 0ES1M, , , , , // V,N .CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 

Case Control Deck: Select desired elements via STRESS^request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose stress is to be 
measured in material coordinates, 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QUADl, (^AD2, TRIAl and TRIA2 elements will appear in the 0ES1M 
output. 

2. A default basic material coordinate system is not provided. For 
this case, set up a rectangular system coincident with the basic 
system. 
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D2. To produce strain/curvature output in both elejnent and material 
coordinates : 

Executive Control Dock: 

1. DIAG 23 

2. ALTER 121 

CURV 0ES1,MPT.CSTO,EST,, / 0ES1M, / C,Y,0UT0PT-O / C,N,1 $ 

0FP 0ES1M,,.,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 

Case Control Deck; Select desired elements via STRESS=request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose strain/curvature is 
to be measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QUADl, QUAD2, TOIAl and ’nUA2 elements will appear in the 0ES1M 
output. 

2. A default basic material coordinate system is not provided. For 
this case, set up a rectangular system coincident with the basic 
system. 
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D3. To produce both strain/curvaturo and stress output in both element 
and material coordinates; 

Executive Control Deck: 

ALTER 121 

CURV 0ESl,MPT,CS'm,EST,, / 0ES1M» / C,Y,0UT0 PThO / C,N,1 $ 

0PP 0ES1M,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C,N»SSST / C,N,+23 $ TURN DIAG 23 0N 

SDR2 CASECC,CSTM,MPT,DIT,EQEX,M,SIL,GPTT,EDT,BGPDT»,,UGV,EST, / 

,,,0ES1A,, / C,N, STATICS $ 

0FP 0ES1A,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

CURV 0ES1A,MPT,CSTM,EST,, / 0ES1AM, / C,Y,0UT0PT«O / C,N,1 $ 

0FP 0ES1AM,,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C,N,SSST / C,N,-23 $ TURN DIAG 23 0FF 

ENDALTER 

Case Control Deck: Select desired elements via STRESS=request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose strain/curvature (or 
stress) is to be measured in material coordinates. 

PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 


2 . 


Remarks : 


1. The selected elements will all appear in the 0ES1 and 0ES1A output. 

2. Only QUADl, QUAD2, 'miAl and TRIA2 elements will generate strain/ 
curvature values in the 0ES1A output. All other elements will 
generate the usual stress values in both sets of output. 

3. Only QUADl, QUAD2, TRIAl and TRIA2 elements will appear in the 0ESIM 
and 0ES1AM output. 
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El. To produce stress output at both elements and grid points (measured 
in material coordinates); 

Executive Control Deck: 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,. // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST,SIL,GPL / 0ES1M,0ES1G / C,Y,0UT0PT»O $ 

0FP 0ES1M,0ES1G,,,, // V,N,CARDN0 $ 

SAVE CARDN0 < 

ENDALTER 

Case Control Deck: Select desired elements via STRESS»request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose stress is to bo 
measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QUADl, QUAD2, TRIAl and TRIA2 elements and grid points connecting 
them will appear in the stress output. 

2. A default basic material coordinate system is not provided. For this 
case, set up a rectangular system coincident with the basic system. 

3. A working core of approximately 3NjN| words (where and N2 are the 
number of cells per side of a rectangular mesh) must be available to 
carry out the interpolation to the grid points of such a mesh. Sub- 
division of the surface mesh into sub-meshes may be necessary for large problems. 
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E2. To produce strain/curvature output at both elements and grid points 
(measured in material coordinates); 


Executive Control Deck: 


1 . 


2 . 


DIAG 23 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,, // V,N,CARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTO,EST,SIL,GPL / 0ES1M,0ES1G / C,y,0UT0PT=O $ 
0FP 0ES1M.0ES1G,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

ENDALTER 


Case Control Deck: Select desired elements via STRESS»request, 


Bulk Data Deck: 

1. The material coordinate system identification number must be coded on 
the MATl or MAT2 card for any element whose strain/curvature is to be 
measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMAP description for details. 

Remarks : 

1. Only QLIADI, QUAD2, TRIAl and TRIA2 elements and grid points connecting 
them will appear in the strain/ curvature output. 

2. A default basic material coordinate system is provided. For 
this case, set up a rectangular system coincident with the basic 
system. 
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3. A working core of approximately (where and N 2 ate the 

number of cells per side of a rectangular mesh) must be available to 
carry out the interpolation to the grid points of such a mesh. Sub- 
division of the surface mesh into sub-meshes may be necessary for 
large problems. 
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E3. To produce both strain/ curvature and stress output at both elements 
and grid points (measured in material coordinates); 

Executive Control Deck: 

ALTER 120,120 

0FP 0UGV1,0PG1,0QG1,0EF1,, // V,N,GARDN0 $ 

ALTER 121 

CURV 0ES1,MPT,CSTM,EST,SIL,GPL / 0ES1M,0ES1G / C,Y,0UT0PT»O $ 

0FP 0ES1M.0ES1G,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C.N.SSST / C,N,+23 $ TURN DIAG 23 0N 

SDR2 CASECC,CSTM,MPT,DIT,EQEXIN,SIL,GPTT,EDT,BGPDT,,,UGV,EST, / 

,,,0ES1A,, / C,N, STATICS $ 

CURV 0ES1A,MPT,CSTO,EST,SIL,GPL / 0ES1A^^0ES1AG / C,Y,0UT0PT=O $ 

0FP 0ES1AM,0ES1AG,,,, // V,N,CARDN0 $ 

SAVE CARDN0 $ 

PARAM // C.N.SSST / C.N.-23 $ TURN DIAG 23 0FF 
ENDALTER 

Case Control Deck; Select desired elements via STRESS=request. 

Bulk Data Deck: 

1. The material coordinate system identification number must be coded 
on the MATl or MAT2 card for any element whose strain/ curvature and 
stress is to be measured in material coordinates. 

2. PARAM 0UT0PT (optional) may be used to select 0FP output device. See 
module DMA? description for details. 
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Remarks : 


1. Only QUADl, QUAD2, TRIAl and TRIA2 elements and grid points connect- 
ing them will appear in the strain /curvature and stress outputs. 

2. A default basic material coordinate system is not provided. For 
this case« set up a rectangular system coincident with the basic 
system. 

3. A working core of approximately 3 NjN| words (where Nj and N 2 are 

the number of cells per side of a rectangular mesh) must be available 
to carry out the interpolation to the grid points of such a mesh. 
Subdivision of the surface mesh into sub-meshes may be necessary 
for large problems. 

Interface with External Post-Processing Programs 

The strain/ curvature (or stress) data generated by NASTRAN as a result of 
this new capability may be picked up by an external program in several ways. 
The simplest is to let j8FP generate punch card images via the case control 
request STRESS (PUNCH) ° . . . 

These card images may then be read by F0RTRAN formatted READ statements in 
the usual way. 

An alternate way is to utilize the utility module 0UTPUT2 to create a 
F0RTRAN binary file which can be directly read by the post-processor program. 
Usage of the utility module 0UTPUT2 and the format of the F0RTRAN-readable 
files are contained in the standard NASTRAN documentation. 
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